Interacting dark energy: generic cosmological evolution for two scalar fields 
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We study the cosmological evolution of two coupled scalar fields with an arbitrary interaction 
term Vr (</>,</?) in the presence of a barotropic fluid, which can be matter or radiation. The force 
between the barotropic fluid and the scalar fields is only gravitational. We show that the dynamics 
is completely determine by only three parameters A^, i — 1,2,3. We determine all critical points 
and study their stability. We find six different attractor solutions depending on the values of Ai and 
we calculate the relevant cosmological parameters. We discuss the possibility of having one of the 
scalar fields as of dark energy while the other could be a scalar field redshifting as matter. 
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A dark energy component is probably responsible for the present stage of acceleration of our universe [1], [2]. Per- 
haps the most appealing candidate for dark energy is that of a scalar field, quintessence [3], which can be either a 
fundamental particle or a composite particle [4] . Within the context of field theory and particle physics it is appealing 
r*"** ' to interpret the dark energy as some kind of particles that interact with the particles of the standard model very 
weakly [5]. The weakness of the interaction is required since dark energy particles have not been produced in the 
accelerator and because the dark energy has not decayed into lighter (e.g. massless) fields such as the photon. It 
is common to assume the interaction between the dark energy and all other particles to be via gravity only, how- 
^v^j ' ever recently interacting dark energy models have been proposed [8]-[ll] . This has been motivated in part by the 
cosmological observations where an equation of state of dark energy may be smaller than minus [1],[2]. In general 
£ — . . fluids with w < — 1 give many theoretically problems such as stability issues or wrong kinetic terms as phantom fields 
[13]. However interacting dark energy [8]-[ll],[12], where the dark energy interacts not only gravitationally with other 
fluids, is a very simple and attractive option which may lead an apparent equation of state smaller than -1 [11]. These 

O fluids can be dark matter, neutrinos or other scalar fields. 
We study in this letter the cosmological evolution of the two scalar fields with arbitrary potentials in the presence 
• of a barotropic fluid, which can be matter or radiation. We show that all models dependence lies on three parameters 
Ph| Xi,i = 1,2,3 defined in eq.(27). We determine the dynamical equations and obtain the attractor solutions as a 
function of these A^ . We find six different attractor solutions depending on the relative size of Ai and we calculate the 
' relevant cosmological parameters. 

This letter is organized as follows. In section II we set up the framework for the cosmological evolution of two scalar 
fields with an arbitrary potential in the presence of a barotropic fluid. In section III we present the conditions on 
the scalar potential for dark energy. In section IV we derive the dynamical first order differential equations and we 
show that the system is determined by only three parameters. In section V we calculate the critical points and we 
study the stability of each solution and we give a few examples. In section VI we study different asymptotic limits 
and we present a discussion on specific particle physics motivated examples. Finally in section VII we present our 
conclusions. 



II. COUPLED SCALAR FIELDS 



Our starting point is a universe filled with two scalar fields </>, ip and a barotropic energy density pb, which can be 
either matter wt, — or radiation uib = 1/3. We will assume that the scalar fields interact via a potential Vt(4>, V 3 ) 
while there is only gravitational interaction between these fields and the barotropic fluid. This work generalizes that 
of a single scalar field and a barotropic fluid [6] . 

One of this scalar fields, namely 0, may be considered as dark energy (quintessence) but it is not necessary to 
interpret <p as dark energy and we will work in a completely general framework. We take the following Lagrangian 
for the scalar fields <j> and ip 

L = ±d^cf>+±d,pd^-V T (0,p) (1) 

with the total potential Vt(4>, </?)■ Instead of working with the total potential Vt we find it useful to separate the 
contribution of the interacting and non interacting terms to differentiate the contribution of the two scalar fields. So, 
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without loss of generality, we take the total potential as 

V T (<t>,<p) = V(<l>)+B(<i>, ( p) (2) 

with V{4>) a potential only depending on (j> and the interacting potential B(<f>, tp) is a function of both fields. Of 
course we could take B(<f>, tp) = h(<f>, tp) + K(tp) with Vr = V(<p) + h(<p, <p) + K(tp) to have a more symmetric potential 
between the fields <j> and tp. However, it is more convenient to keep only the two potentials V and B. 
The equation of motion of </> and tp for a spatially flat Friedman-Robertson- Walker (FRW) universe are 

4> + 3Hj> = -V T4> = -V+ - B (3) 
<p + 3H<p = -V Tv = -B v (4) 

where the subindcx in V and B is defined as V<f, = dV/dcfijB^ = dB/d<f> and B v = dB/dtp. The Hubble parameter 
H = a/a is 

3£f 2 = p = P4, + p v + pb (5) 

where we have taken 8irG = 1 and p is the total energy density, pb the barotropic fluid and p<f>,p v arc defined in 
eqs.(8) and (9). The mass of the scalar fields is given by 

2 d 2 V T 

m = = V<t,<t, + (6) 



d(j) 2 
d 2 V T 
dp 2 



m2 = ^J= B w (7) 



We can also define the energy density and pressure for the fields <f> as 



P4 > = \^ + V{4>), p^\tf -V^) (8) 

and that of tp as 

= 7^P 2 +B{<f>,tp), Pip = ^p 2 - B(<j),p). (9) 
Using eqs.(8) and (9) we can rewrite the dynamical eqs.(3), (4) in terms of the energy densities as 

/?0 + ZHp^l + w<p) = -b = ~B^§ 

Pr + ZHp^l + Wp) = <5 = B (10) 
pb + 3Hp b {l + w b ) = 

where we have included the evolution of the barotropic fluid pb and 

5 = B^^ (11) 

defines the interaction term. The equation of state parameters are given by 

w<b = — = 4-. , Wu, = — = "r — r (12) 

and the time derivative of H is 

H= -^{P4> + Ptp + Pb + P<t> + Ptp + Pb) = ~ (J> 2 + V 2 + Pb(l + Wb)) ■ (13) 



A. Effective Equation of State 



To obtain an effective equation of state we simply rewrite eqs.(10) as 

P4> = ~3H P4> (1 + W ef f), pb = -3Hp b (l + Wbeff) 



(14) 
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with the effective equation of state defined by 

_ _ B^ 

W< p ef f = w 4> + , Wipeff = Wb --^—. (15) 

We see from eqs.(14) that w^ e f f , w ve f j give the complete evolution of p^ and p v . For B^cf) > we have w e ff > w^ 
and the fluid p^ will dilute faster then without the interaction term (i.e. B^cj) = 0) while p v will dilute slower since 
w ipeff > w v . Which fluid dominates at late time will depend on which effective equation of state is smaller. The 
difference in eqs.(15) is [6] 

Aw e ff =W ve ff -Wfeff = Aw-T (16) 

with Aw = w v — Wcj, and T defined as 



T = I P±±P* \ (17) 



3H \ ptppp 
while the sum gives 

^•pWpeff + Q^W^eff = Cl v W v + (18) 

Clearly the relevant quantity to determine the relative growth is given by T and if T > Aw we have Aw e f / < and 
p v will dominate the universe at late times while for T < Aw we have Aw e ff > and will prevail. In the limit 
of no interaction S — B^tfi — we get T = and Aw e ff = Aw > if w$ < w v and p$ will dominate at late times. 
If T = Aw then w ve ff = w 4>eff and the ratio of both fluids p v j p$ will approach a constant value. If the universe is 
dominated by p$ + p v , i.e. 0^, + Q v = 1, then eq.(18) gives 

w<j> < w^eff = w v Q, v + 10^(1 - fi v ) < w v , (19) 

i.e. the effective equation of state is constraint between w<p and w v . 

B. Effective Potential V T 

The effective potential Vt is defined in eq.(2) and we expect the fields (/>,<p to evolve to the minimum of the 
potential, i.e. Vr<j, = V$ + B^ — » and Vr v = B v — » 0. If we want to interpret V as the dark energy potential, 
no fine tuning of the potential V requires that the minimum is at vanishing potential (i.e. V = 0) and V$ < (see 
section III), however including the interaction term B we can now have Vr</> = for a non vanishing total potential 
Vt- A minimum of Vt can be reached if B^ > since V$ < by hypothesis. Taking Vr<p = and the time derivative 
Vt<p = Vt<{><{>4> + Vrjxpip — m 2 <j) + B^^ip — 0, where we used eq.(6), one obtains [11] 

B 4"P l P „ „-3 



with the solution of <p + 3H<j) w giving [14] 



cx a~ 6 (20) 



oca"'*. (21) 



While the derivative of the effective potential Vtc/, is negative the field </> evolves to larger values and even in the limit 
Vt4> — we have with <p > and a positive the interaction term 8 = B^cft. The mass m given by eq.(6) becomes [11] 

m 2 = V M + B H = B H (\ + f * B H ^1 + F^^) (22) 

with T = V^V/Vq (r m > 1 if the field (f> is tracking [3]) and we have approximated p<$, ~ V and p v ~ B in then last 
equality of eq.(22), valid if the kinetic energy is small compared to the potential energy. 
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III. DARK ENERGY 



We may consider the scalar field <j> as the quintessence field, i.e. dark energy. In the absence of any interaction term 
with (p, the lagrangian is simply given by L(<f>) = | <j) 2 — V{4>) and V would be the potential responsible for present 
day acceleration. In this case the slow roll constrains 

1^1 «1, ^«1 (23) 

must be satisfied at present time. As a result of the dynamics, the scalar field will evolve to its minimum and if we 
do not wish to introduce any kind of unnatural constant or fine tuning problem, the minimum of the potential must 
have zero energy, i.e. V\ min = V^\ min — at (f> min [6]. 

In the absence of a interaction B, a finite value <p m in implies that the scalar field (f> oscillates around its vacuum 
expectation value (v.e.v.). If the scalar field has a non zero mass or if the potential V admits a Taylor expansion 
around <fi m in then, using the Hdpital rule, one has limt^oo| V^/T^l = oo and the energy density p^, redshifts with 
u>0 = (n — 2)/(n + 2), i.e. — 0, 1/3 for n = 2, 4 [6]. On the other hand, if (j> m in = 00 then <f> will not oscillate and 
\V^/V\ will approach either zero, a finite constant or infinity. Only in the case |V^/y| going to zero or a constant 
smaller than \f2 will the universe accelerate at late times [6] . 

The absence of an arbitrary dark energy scale requires the slow roll conditions to be satisfied also in the future and 
therefore V is a runaway potential and tends to zero at </> — > 4> m in — oo. From eqs.(23) we see that and also 
approach zero at late times and is therefore negative. Inflation occurs in general for <f> > 1 with a mass m ~ H . 



IV. GENERIC DYNAMICAL ANALYSIS 

To determine the attractor solutions of the differential equations given in eqs.(3) and (4) or (10) it is useful to make 
the following change of variables [6] , [7] 

<t> i fv 

w 1 /b 

and eqs.(10) and (13) become a set of dynamical differential equations of first order 



XlN 



(3+^)*i + y!(A iyI + A 32 , 2 2 ) 



X2N = - I 3 + — I x 2 + W - A 2 y 2 

Vin = --jfyi ~ \J\ Aia;i2/i (26) 
V2N = --rry2 - \ hi (A3X1 + X2X2) y2 



H 2 



H " V 2 
(2x1 + 2x l + ^blb) 



where N is the logarithm of the scale factor a, N = ln(a), 7^ = 1 + w b , f N = df/dN for / = Xi,yi,H (i = 1,2), 
fib = 1 — x\ — x\ — y\ — y\ and 

\i(N) = -^, X 2 (N)^-^, A 3 (7V) = -§. (27) 

Notice that all model dependence in cqs.(26) is through the three quantities Xi(N),i = 1,2,3 and the constant 
parameter -ft, = 1 + Wb- The last eq. of (26) is constraint between — 3 < Hm/H < for all values of Xi,yi and it 
takes the value —3 when the universe is dominated by the kinetic energy x\ + x\ = 1 while it becomes H N /H = 
when the universe is dominate by a constant potential y\ + y\ = 1. The set of equations given in eqs.(26) give the 
evolution of two scalar fields <j), ip with arbitrary potentials in the presence of a barotropic (perfect) fluid with equation 
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of state Wb = 1 — 7b- As mentioned in section II the choice of dividing the total potential Vt into Vr = V{<j>) + B(<fi, tp) 
is without loss of generality and it is convenient in order to distinguish the contribution from both scalar fields. If 
we do not want to consider the contribution from the barotropic fluid we can easily take the limit 75 = in eqs.(26) 
since all contribution form p b is given in H^/H via the term Sl b j b . For Q b ^ we will assume a barotropic fluid with 
< 7& < 2 and 7b = 1 for matter while 7b = 4/3 for radiation. 

We do not assume any equation of state for the scalar fields. This is indeed necessary since one cannot fix the 
equation of state and the potential independently For arbitrary potentials the equation of state for the scalar fields 
w 4> —Pci>/p<p, w v =p v /ptp is determined once 

n* = ^s = + {52=*-* ( 28 ) 

% = ^- 2 =4 + vl ^j = *2-*2 ( 29 ) 

have been obtained. Alternatively we can solve for Xi, y% using eqs. (26) and the quantities w<p = w\ = {x\ — y 2 )/(x 2 +y 2 ) 
and w v = w 2 = (x\ — yl)/( x 2 + 2/1) are > m general, time or scale dependent. In terms of Xi,yt the interaction term 
in eq.(ll) becomes 

8 = = -V63H 3 \ 3 x iy l (30) 
giving an effective equation of state parameters defined in eqs. (15) as 



[2X3^1 x i -vi -y3 A 3 x m ^ 

w *«' = w *-n-ar = n; (31) 

[2 x 3 x iy 2 y / |A 3 x 1 ) 

w *>" = w - + V 3 = % (32) 

and the acceleration of the universe is given by 

- = H 2 + HH N = ~ (fl b (l + 3w b ) + A{x\ + x%) - 2{y\ + y 2 2 )) 

(X A 

= ~(4-3fi 6 (l-«; 6 )-6(y? + y|)) (33) 

where we have used H = HH^ and eq.(26). Clearly acceleration will occur if the universe is dominated by the 
potential y% = (yf + y\) = V T /3H 2 = (V + B)/3H 2 , i.e. for y% > 2/3 - 6 (1 - w b )/2. 



V. CRITICAL SOLUTIONS 



We find the critical solutions to the dynamical equations (26) with x\n = %2N = yiN = V2N = and solve for 
constant values of Aj, i = 1,2, 3. The set of solutions are given in tables I and III. In table I we give, for completeness, 
the unstable critical points and in table II we show the values of the equation of state parameters for these solutions. 
More interesting, we give in table III the attractor (stable) solutions. Constant A implies that the potential V and 
B are exponential potentials, e.g. V ~ e~°"^ with Ai = —V^/V = a constant. However, if the potential is not 
exponential we do not expect to have constant Ai and they will, in general, evolve with time. In this case we can use 
the attractor solution of table III and take the corresponding limit of Aj , which will be either zero, constant or infinity 
[6], to obtain the asymptotic behavior of the solutions. 

In order to determine the stability of the critical points we perturb eqs. (26) around the critical solution Xi,yt and 
we keep linear terms only . The set of eqs. can be written in a matrix form, Zn = MZ, where Z = {8x1,8x2, 8y\,8y2) 
and we diagonalize the matrix M. The stability of the solution requires the real part of all eigenvalues to be negative. 
One of the eigenvalues of model U-I is positive with Ei = 3j b /2 > for 7b > while all other models in table I have 
at least one eigenvalue of the form Ei = 3(2 — g) > for < g < 2. Therefore they are all unstable solutions. 

The eigenvalues of models in table III are given in table IV. We see that depending on the values of Aj the 
eigenvalues can be negative or positive. For any given choice of Ai there is only one stable solution. In table V we give 
the constrains on X's form closure |a;,| < 1, yi < 1 and stability arguments. We see that we have two main conditions 
on Ai. One is the relative size between X 2 and A1A3 while the second is the relative size between \\ + X 2 and A1A3. 
Depending on the relative size of these two conditions the attractor solution will end up either in models S-I,S-II 
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or S-III,S-IV or S-V, S-VI. A further condition on the magnitude of Aj and 37b distinguishes between the different 
models. 

In tables VI and VII we give the values of the relevant cosmological parameters such fi^, Q v , Clb, y\ = y\ + y 2 
and w^jW^jW^ef / ,w ve f / , g = \f2JiXzX\y 2 . As mentioned in section II A the effective equation of state gives the 
correct redhsift of the scalar field when the interaction term is considered. Attractor solutions S-I and S-II are the 
critical solution of a single scalar field in the presence of a barotropic fluid. In these models the energy density p v 
redshifts faster than p§ and pb so w ve ff must be larger than w^eff or Wb- We see from table VII that models S-III 
to S-VI have both scalar fields with the same redshift, i.e. u>0 e // = Wipeff- In models S-III and S-V the redshift is 
equal to that of the barotropic fluid while in model S-IV and S-VI it depends on the value of the different A^ with 
w 4>eff — w <peff < Wb- A universe dominated by the barotropic fluid is possible in models S-I, S-III and S-V, while 
models with no barotropic fluid (Sib = 0) are given by S-II, S-IV and S-VI. An accelerating universe requires y\ to be 
larger than 2/3 — f2&/2 and therefore models S-II, S-IV and S-VI are favored. From table VI and conditions in table 
V we can see that models S-I, S-III and S-V have y\ < 1 — -fb/2 < 1/2 for 7^ > I and therefore do not lead to an 
accelerating universe. 

The eigenvalues of model S-V are 



El V ,l/2 = -j ±7 \l 7772 ( M ) 



Eiv,3/i 




1(2 


- 7fc )(24 7fc 2 [(A 1 - 


-A 3 ) 2 + A 2 ] + (2- 


9 7b )A 2 A 2 ) 


AfA 2 , 


1(2 


-7fc)(87&A 3 [(Ai 


-A 3 ) 2 + A 2 ] + (2- 


-9 76 )AiA 2 ) 



A1A 2 



while of model S-VI are 



AjA 2 , - 3 76 [(Ai - A3) 2 + A 2 ] 

(Aj - A3) 2 + A 2 
A? A 2 - 6[(Ai - A3) 2 + A 2 ] 



ElVIa = 2[(X 1 -X,f + X 2 ] (35) 

p . _ A 2 A 2 -6[(A!-A 3 ) 2 + A 2 ] I ^wl 

Vim ~ 4[(Ai - A3) 2 + A 2 ] ± 4 ((Ai - A3) 2 + A 2 ) 



with 



a v 1 = 2 ABC + 2D 2 - BE - F 2 
A = (Ai - A 3 ) 2 + A 2 , B = X 2 X 2 -6A, C = 3 7fc - A1A3, (36) 

D = 3A(2 + lb ) -2A 2 A 2 , E = A{2X X X 2 - 6(1 + 2 7b )) + 3A 2 A 2 , F = 6A(1 - lb ) + X 2 X 2 

The effective equation of state are given by 

( w 4>eff) = w <P + Jj-, («W/) =W v --£- (37) 

with g given in table VII. For model S-VI we have 

X 2 X\ (6[(Ai - A3) 2 + A 2 ] - A 2 A 2 )(A 2 + A 2 - A1A3) 



(«V) 



VI 



X 2 X 4 2 + (6[(Ai - A3) 2 + A 2 ] - A 2 A 2 )(A 2 + A 2 - A1A3) 



, , A 1 A 2 (2A 1 -A 3 )-6[(A 1 -A 3 ) 2 + A 2 ] 

KV/ (6[(Ai-A 3 ) 2 + A 2 ]-AiA 2 A 3 ) ' (38) 



A. Examples 

We now present four different attractor solution depending on the values of Ai, A2, A3. We show in figures 1-4 the 
evolution of = 0,1,0,^ = f^jfifc and w\ = w ( p,w ( p e ff,W2 = w v ,w ve ff for the different choices of X's. We also 
show the phase space of (27,2/1) and (2:2,2/2) for each case. Since the phase space depends on four variables, namely 
(xi, 2/1,^2,1/2) ! h is no surprising that the curves in the two dimensional space (27,2/1) and (212,2/2) niay cross. 

In figs.l we have A 2 — 5, A 2 = 1, A 3 = 3 and jb = 1 + Wb = 1 which implies that conditions of model S-I in table V 
are satisfied, i.e. 3 = 37h < X\ = 5 and 5 = X\ < A1A3 ~ 6.7. In this case p^ and pb have the same redshift at late 
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Model 


2/i 


2/2 


XI 


£2 


ill = Q$ 


O2 = n v 


fif, 


T T T 

U-l 




















1 


U-II 








±1 





1 








U-III 








4-./1 - T 2 
2 




1 - xi 


xl 





U-IV 








Ai 


1 


6 


1 6 





U-V 








s/6A 3 +A2-v/ A 2+ A 3- 6 


v / 6A 2 -A 3 v' A 2+ A 3- 6 


x\ 


x\ 





Ai^ +A^ 


A 2 + A^ 


U-VI 








^6A 3 -A 2 ^/A^+Ajj-6 


v / 6A 2 + A 3 ^/A^+A^-6 


x\ 







A 2+ A 3 






x 2 



TABLE I: Unstable critical solutions. 



Model 


W<f, 






W ve ff 


T 


U-I 












U-II 


1 




1 






U-III 


1 


1 


1 


1 





U-IV 


1 


1 


1 


1 





U-V 


1 


1 


1 


1 





U-VI 


1 


1 


1 


1 






TABLE II: Equation of state parameters of the unstable solutions of table I. 

times and Sl^/flb approaches a constant value while p v rcdshifts faster with an effective w ve ff > Wb = 0, even though 
u>2 = w v — > — 1, and therefore f2 v tends to zero. The attractor solution has (2:1,2/1) = {y/3/10, \/3/10) ~ (0.55,0.55) 
and (x2, 2/2) = (0, 0) and f2i = 0.6, f2 2 = with fij, = 0.4 and w b — — w ve ff — while w v = —1 and w ve ff = 0.34. 
Since 2/^ = 2/2+2/2— 0.3 is smaller than 2/3 — fi&/2 = 0.87 then from eq.(33) we conclude that there is no late time 
acceleration. 

In figs. 2 we take Ai = 2,A2 = 1/2, A3 = 1/2 and 7$, = 1 + w b = 1 which implies that conditions of model S-IV 
in table V are met, i.e. 0.5 = A 2 + A§ < A1A3 = 1 and 3 = 37b > A2 + A§ = 0.5. In this case p<f, and p v have the 
same redshift at late times and 0,^/0,^ approaches a constant value while fib goes to zero. The effective equation of 
state w<t, e ff and w ve ff have the same late time value w^e// = w ve ff ~ —0.83 while = 1 and w v — —0.91. The 
attractor solution has (2:1,2/1) = (0.2,0) and (£2,2/2) = (0.2,0.96) and Oi = 0.04, ^2 = 0.96 and Of, = 0. In this case 
Vt = 2/1 + 2/i — 0-92 is larger than 2/3 — fib/2 = 2/3 and from eq.(33) we conclude that the universe accelerates at 
late times. 

In figs. 3 we have Ai = 3, A2 = 2, A3 = 3/2 and 76 = 1 + Wb = 1 which implies that conditions of model S-V in table 
V are satisfied, i.e. 6.25 = X\ + A 2 > X 1 X 3 = 4.5, 9 = A 2 > A1A3 = 4.5, 18.75 = 37&[(Ai - A3) 2 + A 2 .] < A 2 A| = 36 and 
9 = A 2 > 37 2 /2 = 3/2. In this case p^, p v and p b have all the same redshift at late times and fi^, fi v , fi b approach a 
constant value. The effective equation of state w^, e ff and w ve ff have the same late time value w^, e ff — w ve ff — w b — 
while — 0.4 and w v = —1/3. The attractor solution has (2:1,2/1) = (0.4,0.27) and (#2,2/2) = (0.3,0.43) and 
fii = fi = 0.24, fi 2 = fi v = 0.29 and fi b = 0.47. Since y\ = y\ + y\ = 0.26 is smaller than 2/3 - fi 6 /2 ~ 0.43 then 
from eq.(33) we conclude that the universe does not accelerate at late times. 

Finally, in figs. 4 we take Ai = 1, A 2 = 1, A3 = —3 and jb — 1 + Wb = 1 which implies that conditions of model S-VI 
in table V are satisfied, i.e. 10 = A^ + A§ > A1A3 = -3, 1 = X\ > A1A3 = -3, 51 = 37 6 [(Ai - A3) 2 + A 2 ,] > A 2 A 2 . = 1 and 
17\/6 = V6|Ai|[(Ai - A3) 2 + A 2 ] > A 2 A 2 = 1. In this case p</> and p v have the same redshift at late times and fi^/fi^, 
approaches a constant value while fib goes to zero. The effective equation of state w^ e ff and w ve ff have the same late 
time value w^eff — w<p e ff — —0.98 while w,f, = —0.99 and = —0.92. Notice that at some redshifts the effective 
equation of state take values w^ e ff < —1. The attractor solution has (2:1,2/1) — (0.02,0.87) and (2:2,2/2) = (0.1,0.48) 
and fii = fi = 0.76, fi 2 = % = 0.24 and fi fc = 0. Since y\ = y\ + y\ = 0.99 is smaller than 2/3 - fi fc /2 = 2/3 then 
from eq.(33) we conclude that the universe does accelerate at late times. 
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Model 


2/i 


y-2 


Xl 


X2 


S-I 

S-II 
S-III 

Q TAT" 
b-1 V 

S-V 
S-VI 


\ 


/3(2- 7i ,) 7 (, 

V 2 ^ 

V 6 




/3(2- 7b )7i,(A2 + A2-A 1 A 3 ) 





/3(2- 7b )76 

V 2 (^+^) 

V I 6 

/3(2- 7b )7 6 (A 1 -A 3 ) 


V 2 Ai 
Ai 

/3" 7b A 3 

A 3 
76 

AlA 2 , 





/3" 76 A 2 

A 2 
76 

/T7b(Ai-A 3 ) 


/ 2A 2 A 2 


V 2A 1 A 2 


V 2 AjA 2 
A 1 A 2 (A 1 -A 3 ) 




-A 3 ) 2 + A^]-A^A^))(A^ + A^-A 1 A3) 


^/A^As-AOtA^A^-SKAi-Agja+A 2 ]) 


VB((Ai-A 3 ) 2 +A|) 


s/6((Ai-A 3 ) 2 +A 2 ) 




n/6((Ai — A 3 ) 2 + A 2 ) 



TABLE III: Stable critical solutions. 



Models 
S-I 
S-II 
S-III 

S-IV 
S-V 
S-VI 



Eh 



1(70-2) 

5(6 — A 2 ) 

1(7,-2) 

i(A| + A|-6) 
£iv,i 



£i 2 



§76(1 "XT 



1(6 -A?) 

3 7 b(A 2 ' + A 2 -AiA 3 ) 



I(A|+A|-6) 

S«V,2 

Eiyi,2 



1(70-2) + - 



1(76 



2) + | 



(2-7 b )(24 7 2 + (2-976)A 2 ) 



-3 76 + Ai 



(2- 7b )(24 7 2 +(2-97 b )(A 2 + A 2 )) 
T 



A 2 +A 

A2 + A| - 37b 
Eiv,3 
Eivi,z 



[(70-2) 



!(7o-2) 



gi 4 

(2-76)(24 7 2 + (2-976)A 2 ) 



^Ai(Ai — A3) 



i(A| + A3 — A1A3) 

Eiv,4 
EiviA 



TABLE IV: Eigenvalues for the different stable solutions of table III. Stability requires the real part of all eigenvalues to be 
negative. Eigenvalues Eiv are given in eqs.(34) while eigenvalues Eivi in eqs.(35). 



Model 


Constraint 


Constraint 


S-I 
S-II 
S-III 

S-IV 
S-V 

S-VI 


A 2 > 3 7b 
A 2 < 3 76 
A 2 + A 2 < A1A3 
A 2 + A3 < A1A3 
A| + A§ > A1A3 

A? > 3 7fc 2 /2 
A| + Ai > A1A3 
\/6|Ai|[(Ai - A 3 ) 2 + A|] > A?Ai 


Ai < A1A3 
A 2 < A1A3 
\l + A| > 3 76 
\l + A| < 3 7 6 
A 2 > A1A3 
3 7 o [(Ai - A 3 ) 2 + A|.] < A 2 A 2 

A 2 > A1A3 
37o[(Ai-A 3 ) 2 + A|] > A 2 A| 



TABLE V: Constrains on Ai, A2, A3 for the different stable solutions of table III from closure and stability arguments. 



Model 




Q2 = 


fib 


= y? + !/i 


S-I 

S-II 
S-III 

S-IV 
S-V 

S-VI 


1 

■ill A 3 





3 7 6(2A 2 + (2- 7b )A 2 ) 


1 376 



^ 376 


3(2-76)76 
2A 2 

1_£ 
x 6 

3(2-76)76 


2(A 2 +A 2 ) 2 
\ 2 
±3. 
6 

37b(2A 2 -A 3 (2- 7b )(A 1 -A 3 )) 


2(A 2 + A 2 ) 2 

\ 2 
1 A 3 

1_ 6 

37b(Ai-A 3 )(2A 1 - 7b A 3 ) 


A 2 +A 3 



n 3 7b ((A 1 -A 3 ) 2 + A 2 ) 


2(A 2 +A 2 ) 

■1 A 2 + A 3 

6 

3(2-7b)7b[(Ai-A 3 ) 2 + A 2 ] 


2A|A 2 

(6[(A 1 -A 3 ) 2 +A 2 ]-A 2 A 2 ))(A 2 + A 2 -A 1 A 3 ) + A 2 A| 


2A 2 A 2 

Al(A 1 -A 3 )(6[(A 1 -A 3 ) 2 +A 2 ]-A 1 A 2 A 3 ) 


1 A 2 A 2 
A 1 A 2 




2A 2 A 2 
A A 2 

1 A 1 A 2 


6[(A!-A 3 ) 2 + A 2 ] 2 


6[(A 1 -A 3 ) 2 +A 2 ] 


6[(A 1 -A 3 ) 2 + A 2 ] 



TABLE VI: Energy density for the different stable solutions of table III. 
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Model 
S-I 
S-II 
S-III 

S-IV 
S-V 
S-VI 



-1 + 4f 



1 4- 2lbX 2 

^ 2A a , + (2-7 b )(A 3 -A 1 )A 3 



-1 + 



2A|+(2- 7i) )A| 
-1+ 2A ° 



6" A 3 
2(l-7b) A i+7b A 3 
7b X 3— 2X 1 

(w v )vi 



-1 



W <t>eff 
W b 

A 2 

x T 3 

aS+a 2 , 



"'6 



1 J A 1 A 2 

_1 3[(A! — A 3 ) 2 + A^] 



w <peff 

> Wb,W<t,eff 

> W b ,W^ e ff 

A2 + A? 



-1 + 



1 I A 1 A 2 

_1 SKAi-Aap+A?] 



2 ( A I+ A f) 
A 2 (A 2 + A 2_ 6) 



3(2-76)7^3^! -A 3 ) 
2A^ 

AjA|A3(Ai-A3)(6[(Ai-A3) 2 + A|]-Afo|) 



TABLE VII: Equation of state parameters of the stable solutions of table III. The equation of state parameters wvi are given 
in eqs.(38). 




2 4 



10 12 14 



. 8 



0.6- / 



0.4 



0.2 




-0.5 





Xl 



0.5 



0.5 



-0.5 




2 4 6 8 10 12 14 
N 




-1 -0.5 0.5 

X2 



FIG. 1: We show the evolution of Q,f, = = 02,Ob (blue (solid), red (dotted) and black (dashed), respectively) and the 

equation of state parameters w^, = wijW^eff (blue (solid), blue (dotted)) and w v = u>2,w ve ff (red (solid), red (dotted)) as a 
function of N = Log[a], for A? = 5, A2 = 1, A3 = 3 and 7(, = 1 + Wb = 1. With these choice of A's the attractor solution is given 
by model S-I. The attractor solution has (21,2/1) = (y/3/10, a/3/10) ~ (0.55,0.55) and (£2,2/2) = (0,0). 



VI. ASYMPTOTIC BEHAVIOR 



Special cases can be studied by taking different limits of the parameters A. From tables V and VII we can determine 
which values of A^ are required for any particular late time behavior we wish to have. For example if we are interested 
in the behavior of having only two scalar fields and no barotropic fluid (£!{, = 0) we can take the limit 7^ = of models 
in table III. Only models S-II, S-IV and S-VI survive and 2^,2/4 take the same values as in table III. If we prefer a 
universe dominated by the barotropic fluid than models S-I, S-III and S-V need to be considered. An accelerating 
universe requires y\ to be larger than 2/3 — fib/2 and therefore models S-II, S-IV and S-VI are favored. 

If we take the limit Ai = then only models S-II and S-VI remain consistent and they have (xi, yi) = (0, 1) and 
(#2)2/2) = (0,0). Models S-I and S-V do not satisfy the closure condition \xi\ < l\,yi < 1) while models S-III and 
S-IV are no longer stable if A 2 / / A 3 (c.f. eigenvalues of table IV). In the limit Ai — * 0, the first derivative of the 
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-1 -0.5 0.5 1 -1 -0.5 0.5 1 

Xl x 2 



FIG. 2: We show the evolution of = Qi,fl v — Q2,Qb (blue (solid), red (dotted) and black (dashed), respectively) and the 
equation of state parameters w^, = wi,w$ e ff (blue (solid), blue (dotted)) and w v = W2,'W V efS ( re d (solid), red (dotted)) as a 
function of TV = Log[a], for Ai = 2, A2 = 1/2, A3 = 1/2 and ji, = 1 -+- Wb = 1. With these choice of A's the attractor solution is 
given by model S-IV. The attractor solution has (21,3/1) = (0.2,0) and (22,1/2) = (0.2,0.96). 



Model 


yi 


2/2 


2l 


X2 


S-I 


/3(2— no-y,,) 

V '^i 





/Tib 
V 2 Aj 





S-II 
S-III 


\/»-* 






/3(2- 7i ,) 7b ) 
V 2A 2 







/3~7b 
V 2 A 2 


S-IV 










A 2 


S-V 


/3(2- 7b ) 7b 

V 2A i 


/3(2- 76 ) 7i , 
V 2A 2 


/"3 7b 
V 2Aj 


/3" 7 b 
V 2 A 2 


S-VI 




^/A2( 6 [(A2+A2]-A2a2) 




A?A 2 




v'efAf + A^) 


v^A^ + A^) 





TABLE VIII: Attractor solutions of table III in the limit A 3 = 



potential approaches zero faster than the potential itself and examples of this kind of behavior are given by potentials 
of the form V = V (j)~ n ,n > 0. 

For A2 = only model S-V is no longer consistent. Models S-I and S-II remain the same as in table III and model 
S-III becomes (xi, V l ) = ( ^3 /2 7b /A 3 , 0) and (x 2 ,y 2 ) = (0, (^3(2 - 7b )V2/A 3 ), S-IV is now {x uVl ) = (A 3 A/6,0 ) 

and (£2,2/2) = (0, (y/1 - A§/6) while S-VI becomes (x u Vl ) = (0, y/X 3 /(\ 3 - Ai)) and (a*, y 2 ) = (0, ^i/(Ai - A 3 )). 

In the case A 3 = all models remain valid from closure arguments and are given in table VIII. However, models S-I 
to S-IV are no longer stable for Ai 7^ 7^ A 2 , as seen from table IV, and only models S-V and S-VI are the attractor 
solutions. 

Finally, if we impose the condition Xiyf + A 3 y| = 0, which minimize the potential as a function of <j>, i.e dVr/dcj) = 
d(V + B)/d(j) = 0, there are only three attractor solutions to eqs.(26). The first case is model S-II {x\,y\) — (0, 1) 
and (0:2,2/2) = (0,0) (c.f. table III in the limit Ai = — A3J/I/2/1 = 0)- The other two models are S-III and S-IV of table 
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0.5 



. 5 
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0.2 
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\ \ 


/ / / 




' ' / / / ^ 


J X 


1/ / / jf 
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\ / y y /(V 
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\ 

1 



-0.5 





x 2 



0.5 



FIG. 3: We show the evolution of Q^, — f2i,f2 v = Q2,Qb (blue (solid), red (dotted) and black (dashed), respectively) and the 
equation of state parameters w$ = w\,w$ e ff (blue (solid), blue (dotted)) and w v = W2,w ipe {{ (red (solid), red (dotted)) as a 
function of N = Log[a], for Ai = 3, A2 = 2, A3 = 3/2 and ji, = 1 -+- Wb = l.With these choice of A's the attractor solution is 
given by model S-V. The attractor solution has (air, 2/1) = (0.4,0.3) and (22,1/2) = (0.3,0.4). 

Ill with the limit A3 = —\\y\jy\ — (c.f. table VIII). There is no stable solution with Ai, i — 1,2,3 constant and 
2/12/2 7^ 0. Let us take the limit \\y\ + A3J/5 = in eqs.(26). The evolution for x\ is 



Xl 



(3 + ^f)<0, 



(39) 



and since —3 < H^/H < for all values of Xi, yi and we conclude that x\ will approach its minimum value (i.e. 
x\ — > 0). If |Ai| < 00, so that A3X1 — > 0, then the evolution of y\ becomes 



yiN 
2/i 



Hn 



Ai x x 



> 



(40) 



and ?/i will increase to its maximum value (i.e. y\ — > 1). Since x\ + y\ + x\ + y\ + fit, = 1 then for yx = 1 we get 
xi = X2 = 2/2 = ^6 = giving model S-II. 



A. Particle Physics Model 



Let us now take a specific choice of potentials motivated by particle physics. We consider a factorisable interaction 
potential 



B(<f>,v)=h(<l>)F(<p). 

In this case the Ai parameters become only functions of a single field 

Atw=-m A 2( ,)=-m w)- m 



v(4>) 



K4>) 



(41) 



(42) 
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FIG. 4: We show the evolution of 0^ = — Q2,Qb (blue (solid), red (dotted) and black (dashed), respectively) and the 

equation of state parameters w^, = wi,w$ e ff (blue (solid), blue (dotted)) and w v = W2,'W V efS ( re d (solid), red (dotted)) as a 
function of N = Log[a], for Ai = l,Aa = 1, A3 = —3 and ji, = 1 -+- Wb = 1. With these choice of A's the attractor solution is 
given by model S-VI. The attractor solution has (xi,yi) = (0.02,0.87) and (2:2,2/2) = (0.1,0.52) 



where the prime denotes derivative w.r.t. its argument. Let us take a simple example where the potential for the 
scalar field <fi has an exponential behavior, widely used in particle physics as for dark energy potential, while we take 
a power law potential for ip which would represent a standard scalar field. It is well known that for a single scalar 
field with potential of the form <p n the energy density redshifts as w v = (n — 2)/(n + 2) giving w v = 0, 1/3 for n = 2, 4 
[6]. So, we take V = Voe'^.B = B e^<p n with V a > 0, B Q > 0. In this case Ai = a, A 3 = -/?, A 2 = -n/ip . The 
effective potential defined in eq.(2) is minimized for 

V T(j> = V (/> + B^ = -aV + 0B = O (43) 
V Tv = B v = nBoe-^cp 71 - 1 = (44) 

which implies V/h = {V a / B )e- (a+l3 ^ = (f3/a)tp n , i.e. 

= -Log[A V n }/(a + (3), h = e -W(«+moB[A*p«] = A -p/( a +e) ^ (45) 

with A = B j3/V a and A should be positive and we take a > 0, (3 > 0. The condition hip 71 ' 1 — ► becomes 

V" 1 = /^-/3/(«+/3y(«-l)«+«/3]/(«+/3) (46) 

and for \(n — l)a + n0\/(a + 0) > (valid for n > 1) we find ip — > giving |Aa| — > 00. From table V we see 
that for I A2 1 S> I Ai|, I A3 1 the stable solution are given by model S-V or S-VI depending whether 37^ is smaller or 
larger than \\ = a 2 . Model S-V has (fi x , 2 , Vt) = (37b/A?,0, 1 - 37 h /A?,3(2 - 7 h )7 fc /2A?) and model S-VI has 
(Oi, SI2 , j/y) = (1, 0, 0, 1 — Af/6). We see that in both cases the amount of energy density from the field ip vanishes 
and the solutions reduce to a single scalar field <j> which depending on the value of Ai = a we can have a universe 
completely dominated by the scalar field (f> with an accelerating universe for Aj < 2 or a scalar field redshifting as the 
barotropic fluid for \\ > 3j b - The interaction term g given in table VII vanishes for models S-V, S-VI in the limit 
I As] — > 00. The equation of state for <p for model S-V is w v = (2(1 — 7&)Ai + 7bA 3 )/(7bA3 — 2Ai) while for model S-VI 
we have w v = (6 + A1A3 - 2A?)/(AiA 3 - 6) (c.f. eq.(38)). 
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VII. SUMMARY AND CONCLUSIONS 

We have studied the dynamical system of two scalar fields with arbitrary potentials in the presence of a barotropic 
fluid in a FRW metric. We have shown that all model dependence is given in terms of three parameters, namely 
Xi(N) — — V^/V, A 2 (iV) = —B v /B, A 3 (A r ) = —B^/B, and we have calculated all critical points and determine their 
stability. We have seen that there are six different attractor solutions given in table III. For a given choice of \ the 
attractor solution depends on the relative magnitude of Ai , A2 , A3 and 7^ as discussed in table V. We have calculated 
the relevant cosmological parameters and we have shown the phase space for four different attractor solutions. Finally 
we have discussed the different asymptotic limits of the parameters Ai and we have studied a special type of scalar 
potential motivated by particle physics. 
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